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We argue that the shape of the dispersion along the nodal and antinodal directions in the cuprates 
can be understood as a consequence of the interaction of the electrons with collective spin excitations. 
In the normal state, the dispersion displays a crossover at an energy where the decay into spin 
fluctuations becomes relevant. In the superconducting state, the antinodal dispersion is strongly 
affected by the (tt, tt) spin resonance and displays an 5— shape whose magnitude scales with the 
resonance intensity. For nodal fermions, relevant spin excitations do not have resonance behavior, 
rather they are better characterized as a gapped continuum. As a consequence, the S— shape 
becomes a kink, and superconductivity does not affect the dispersion as strongly. Finally, we note 
that optical phonons typically lead to a temperature independent S— shape, in disagreement with 
the observed dispersion. 

PACS numbers: 74.25.Jb, 74.72. Hs, 79.60.Bm, 74.50.+r 



I. INTRODUCTION 

Angle resolved photoemission (ARPES) experiments 
are valuable sources of information about the shape of 
the Fermi surface in the cuprates, and the frequency, mo- 
mentum and temperature dependence of the electron self- 
energy. The subject of this paper is an analysis of the dis- 
persion along various momentum cuts. These dispersions 
have been obtained by several groupsiiSi2i^iSiSiL2i2ii£ by 
high precision measurements of the momentum distribu- 
tion curves (MDCs), by which the spectral function is 
obtained at a given energy by making scans along direc- 
tions normal to the Fermi surface. The spectral function 
v4(k,Li;) — (l/7r)|/mG(k, w)! is related to the self-energy 
Ek('^) as 
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(1) 



Near the Fermi surface, et ~ vp{'kp){k± — kp) where 
vfO^f) is the bare value of the velocity. There are sev- 
eral reasons (both theoretical and experimental) to be- 
lieve that the self-energy weakly depends on the value 
of k± normal to the Fermi surface, and can be approxi- 
mated as w Ekjj (w). For a given kj? specified by a 
cut, the MDC spectral function A{k±,uj — const) is then 
a Lorentzian centered at w_F(kF)(fc_L — kp) — uj — 'Si^p{u>) 
with a HWHM equal to J:'^^{uj)/vp{kF) (Ref.H. In a 
generic Fermi liquid, the self-energy is linear in lo at the 
lowest energies: T,\^^{iu) = —Xif^pCu. The position of the 
MDC peak then determines the renormalized Fermi ve- 
locity Vp = vp/{l + Akjj). At higher energies, A becomes 
frequency dependent, and the dispersion deviates from 
the linear form. 

The MDC data have revealed several characteristic fea- 
tures of the dispersion which need to be explained: 

• In the normal state, the dispersion along both 
the nodal and antinodal directions shows a rela- 
tively smooth crossover from a linear behavior at 



small binding energies to a more steep behavior 
above roughly 50 — 70meV. This effect has been 
observed in Bi221'2ii2Mi^, Bi220lM, Bi2223&, 
Na^CCOC^, and LSCO^. In the last case, the 
crossover is sharper and more resembles a kink. 

• The renormalized Fermi velocity Vp along the nodal 
direction weakly depends on doping^ and in _Bz2212 
equals l.QeVA^. At the same time, at high ener- 
gies, the dispersion is strongly doping dependent, 
becoming steeper with underdoping^. 

• In the superconducting state, the dispersion along 
the antinodal direction develops an 5*— shapeii, 
with a "negative" dispersion between 60 and 
SOmeV. This feature develops with decreasing tem- 
perature in an order parameter like fashion, with an 
onset temperature at Tc for overdoped samples, and 
somewhat above Tc for underdoped samples 
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• The dispersion along the nodal direction does not 
develop an S*— shape below Tc- Instead, in the 
superconducting state the crossover gets sharper, 
with a kink-like feature near 70meV developing in 
gii2212^i3,4,5,8 aj^j Bi2223S,. This extra "sharp- 
ness" has a temperature dependence similar to that 
of the antinodal dispersion mentioned above^. In 
Bi2201 and LSCO, on the other hand, the nodal 
dispersion does not change much between the nor- 
mal and superconducting states^. 

• The high energy nodal dispersion never recovers to 
the bare dispersion. It remains linear to the highest 
binding energy studied^, with an interpolation to 
a k point at zero energy displaced well inside the 
Fermi surface. 

Theoretical scenarios proposed to explain the data dif- 
fer primarily on whether the electron-electron or the 
electron-phonon interaction is responsible for the ob- 
served behavior. In the electron-electron scenario, the 
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crossover from a linear dispersion at the lowest ener- 
gies to a more steep dispersion at higher energies has 
been identified^ with the crossover from Fermi liquid 
to non-Fermi liquid behavior. In the superconducting 
state, the 5— shape dispersion along the antinodal di- 
rection has been associatedii'i^'i^^ with the interaction 
with the (tt, tt) spin exciton, which below Tc emerges at a 
frequency oJres < 2A due to a feedback of the pairing on 
the spin susceptibility. The interaction with the exciton 
gives rise to a which strongly increases and then 

rapidly drops as lo approaches = A -I- u>res ■ This gives 
rise to the observed 5— shape dispersion. In optimally 
doped Bi2212, both A and cOres are close to AOmeV^, 
i.e., ujQ ^ 80mey. The magnitude of the 5— shape piece 
is stronger for |cli| < ujq, since for \uj\ > uq, S" rapidly 
increases. Both the value of loq and the asymmetry of 
the S*— shape dispersion agree with the dataiSiiMi. 

For nodal fermions, scattering by qo = (tt, tt) shifts 
the nodal Fermi point to an energy about 0.7eV above 
Ep. This energy is too high to expect any apprecia- 
ble effect on the low energy dispersion. Still, Eschrig 
and Normarti^ argued that the kink near 70meV can 
be explained by the interaction with the spin resonance, 
as the resonance has a finite momentum width around 
(tt, 7r)ii. They used a phenomenological form of the spin 
susceptibility with a sharp S function in frequency at ujres 
and a Lorentzian in momentum space with a width of 2 
lattice constants, with the momentum smearing giving 
rise to resonance scattering between the node and other 
Fermi surface points. But a good description of the nodal 
fermion spectrum required the inclusion of a gapped con- 
tinuum in the spin excitation spectrum, which acts to 
smear the S*— shape into a kink, and also gives rise to 
the linear uj behavior in /toE observed at higher binding 
energies. 

The electron-phonon scenario for the dispersion was 
put forward in Ref. 0. The key difference with the 
electron-electron scenario is in the interpretation of the 
normal state ARPES data - Ref. argued that there 
is a sharp kink (rather than a crossover) in the disper- 
sion in the normal state, and that the kink energy is 
about the same in all materials studied {LSCO, Bi2201, 
and Bi2212). They further argued that the kink effect 
is rather isotropic in the Brillouin zone (in disagreement 
with other work^iSiiii) . They speculated that this simi- 
larity implies that superconductivity plays a secondary 
role in the phenomenon, and that the features in the dis- 
persion can be reproduced by coupling an electron to a 
bosonic mode unrelated to superconductivity. They sug- 
gested a (tt, 0) optical phonon with an energy of 55meV— 
as the best candidate. 

In this paper, we distinguish between these two possi- 
bilities and argue in favor of a spin-fluctuation scenario. 
We first analyze the spin-fluctuation scenario in more 
detail. We argue that the spin resonance scattering is 
effective in scattering antinodal fermions near the Fermi 
energy, but is not effective for nodal fermions since the 
bosonic momenta which connect a nodal point with other 



points on the Fermi surface are far removed from (tTjTt). 
Rather, the most effective low-energy scattering for a 
nodal fermion is to the antinode, where the density of 
states has a singularity. This scattering gives rise to a 
kink in the self-energy of a nodal fermion at ^ 2A which 
generates a kink in the dispersion at the same energy. 
We flnd that the resonance scattering emerges away from 
the nodal direction, and the magnitude of the resulting 
5— shape dispersion progressively increases as the antin- 
ode is approached^°. 

We then argue that the interaction with an Einstein 
phonon gives rise to a temperature independent 5— shape 
dispersion for all cuts normal to the Fermi surface. This 
is difficult to reconcile with both the antinodal dispersion, 
for which the S*— shape is present but only emerges below 
Tc, and the nodal dispersion which does not display an 
5— shape form at any temperature. 



II. 



MAGNETIC SCATTERING 



In the magnetic scenario, the fermionic self-energy 
originates from the strong spin component of the 
electron-electron interaction in the particle-hole channel 
and can be viewed as coming from scattering by collec- 
tive spin fluctuations. To lowest order, the corresponding 
self-energy is given by 



Skp(w) = - 
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d^qdnGl.jLu + n)xs{ti.n) (2) 



where gs is the spin-fermion coupling. Here G'^ is the 
bare Green's function (in the normal state, G'°(k, w) = 
1 / {uj — VFO^F){k± — kp))), and Xs{<l, t^) is the dynamical 
spin susceptibility for which one has. 
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where Xs (q) is the static part of the susceptibility which 
is believed to be peaked at or near the antiferromagnetic 
momentum qo = (tt, tt), and nq(ri) (subject to nq(0) = 
0) accounts for the spin dynamics and is proportional to 
the dynamical part of the full particle-hole bubble. For a 
non-diagonal q that connects points on the Fermi surface, 
the polarization operator has the forr^ 



n.(f^) = *EY / 



cLu 



1 - 



A+ A_ -I- uj+uj^ 



-A2 



- A2 



(4) 

Here u!± = oj ± il/2 and A± — A(k,„ ± q/2) are the 
values of the c?— wave gap at the points km ± q/2 which 
are simultaneously on the Fermi surface, the summation 
m being over a discrete set of these points. The prefactor 
7q depends on the coupling gs and the angle between the 
Fermi velocities at k„j ± q/3i^. In principle, the pairing 
gap A depends on frequency, but this dependence is not 
essential for our purposes and we neglect it for clarity. 
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A. Normal state 

1. polarization operator 

In the normal state A = 0, and nq(ri) is purely imag- 
inary: nq(f7) = Z7q|f7|. This is the expected result as 
once q is such that two k— points separated by q can 
be simultaneously put on the Fermi surface, the polar- 
ization bubble contains a Landau d amping term . This 
term generally has the form y^^vpq)-^ — ft^ but in 
our case q is finite and vpq well exceeds fi. The true po- 
larization bubble also contains a frequency independent 
piece, but this piece comes from fermions with energies 
comparable to the bandwidth and is already incorporated 
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into X7^{^)- Note that this separation is consistent with 
a Kramers-Kronig (KK) analysis: a KK transformation 
of /mn(0) = 7^ does not produce a universal piece of 
ReIl{Q) independent of the upper limit of the frequency 
integration. 



2. fermionic self-energy 

Substituting the relaxational Xs(q, into the self- 
energy, introducing a small q via q q + Qo and lin- 
earizing the fermionic dispersion near the Fermi surface, 
we obtain from Q 



J dqi\dq±dn 



u + ci - vfO^f + q + qo)<7_L 
I 



\q + qo) -i7q+qo|f^l 



(5) 



where N stands for normal state. For consistency with 
the assumption that the self-energy weakly depends on 
Ck, we assume that the fermionic propagator changes 
much faster with q± than the bosonic x{q, ^^), i-e., that 
the Fermi velocity is much larger than the "spin" veloc- 
ity. We then integrate over momentum q± normal to the 
Fermi surface only in the fermionic propagator, and set 
qj_ in the bosonic propagator equal to its value at a dis- 
tance between kj? and some other point on the Fermi 
surface which is parametrized by g|| . The integration 
over q± is straightforward, and performing it using the 
fact that Xs(q, il) is an even function of frequency, we 
obtain 
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vpCkp + qo + qii) 



Xs (qo +q||) - i7qo+q||^^ 



(6) 



The remaining integral over depends on the momen- 
tum dispersions in Xs(q) and vpCkp) along the Fermi 
surface, both are inputs for the low-energy theory. For 
a qualitative understanding of the crossover in the dis- 
persion, we assume momentarily that Xs(qo + q||) is flat 
near qo, and that 7qo+q|| and vfO^f + qo + q||) are also 
momentum independent. We then immediately obtain 
from ® 



i\uj\_ 

UJsf 



sgn{uj) (7) 



where A — ig'^Xs/i^TTVp) and ujsf — {Xsl)~^ (we use the 
same notation as in earlier worksiSiSS). Separating real 
and imaginary parts of the complex logarithm, we obtain 



from © 



i?eEjvkF(w) = — Aws/arctan 

/mSAr,k^(tj) = -A^^ log 1 + - 
2 oj 



sgn{uj) (8) 



At low frequencies, one indeed recovers Fermi liquid be- 
havior: 



. (tj) — — A [LiJ + i 



. uj\uj\ 



2uj 



(9) 



On the other hand, at frequencies larger than Wsf, the 
self-energy nearly saturates 

EV(t^) « -{TT/2)Xujsfsgn{cu), S'^(cj) (x log \cu\ (10) 

The evolution of with frequency gives rise to a 

crossover in the normal state dispersion lo — — 
VF{k± — kp) around lo — tOsf- We illustrate this in Fig^ 
We clearly see that the dispersion is linear below ojsf, 
with the eff'ective velocity v*p — Wi?/(1 + A). However, 
above uJsf, the dispersion crosses over to a more steep 
form which also yields an intercept at a finite k± — kp 
if extrapolated formally to zero energy. This crossover 
behavior is consistent with the one observed experimen- 
tally. Note also that is almost linear in frequency 
in a relatively wide frequency range above LUsf. This 
quasi-linearity seems to be a generic property of in 
the crossover region between Fermi liquid behavior at 
small frequencies and quantum-critical, non-Fermi-liquid 
behavior at larger frequencies. 

The inclusion of the momentum dependences of Xs (q) , 
7q and vp(kp) gives rise to some uj dependence of S^(aj) 
at high frequencies, and to the angular dependence of 
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FIG. 1: The self-energy and the dispersion in the normal 
state for a flat static susceptibility Xai'i) near qo = (tt, tt). a) 
ReTiN, b) and c) Ira's n, with c) over intermediate frequencies 
where ImTiN{^) displays a quasi-linear behavior, d) disper- 
sion uj — ReTiiu)) — vpAk (Ak = k — kp), with the dashed 
line the free fermion dispersion. The coupling is A = 2. 



the coupling constant A, but the crossover near ujsj 
still survives. To illustrate this, in Fig[21 we plot the 
the dispersion obtained for the Ornstein-Zernike form 
of xJ^(q) oc 1 + (q^)^ with constant 7 and vp- We 
see that the dispersion is again linear at small w with 
Vp = vp/i^ + Aky), and crosses over to a more steep 
dispersion above cOsfO^p)- Observe also that S"(a;) is 
again nearly Hnear above uisfOf-r)- The crossover fre- 
quency 0Jsfi)<.F) = (7Xs(k_F))~^ is smallest for an antin- 
odal fermion and largest for a nodal fermion simply 
because the node-node distance is smaller than qo — 
(tt, tt), and hence for node- node scattering xJ^(q) ^ 
(1 -I- (|q — qo 1^)^). For an antinodal fermion, on the other 
hand, the antinode-antinode scattering involves momenta 
very close to qo, hence xl^i'i) i^ smaller as the piece 
(|q — lolO^ i^ absent. This effect is, however, partly 
compensated by the fact that 7q is enhanced around a 
nodal point and formally diverges for node-node scat- 
tering because the Landau damping blows up when the 
velocities of the two fermions in the particle-hole bub- 
ble become antiparallel to each other, as is the case for 
nodal fermions^. Previous calculations showSi that, as 
an interplay between the two effects, the variation of cosf 
along the Fermi surface near optimal doping is relatively 
modest, i.e, the crossover frequency for the normal state 
does not vary substantially along the Fermi surface. 



B. Superconducting state 

1. polarization operator 

We begin with the polarization operator, Eq. (01 . Ap- 
plying the spectral representation to one can imme- 



FIG. 2: Same as Fig. Q but for the Ornstein-Zernike form of 
the static spin susceptibility. In both cases, the crossover oc- 
curs around lo — u^f The dashed line in d) is the free fermion 
dispersion. The dashed-dotted line in d) is an extrapolation 
from high frequencies. The extrapolated dispersion crosses 
the vertical axis at a negative Ak, i.e., for k inside the Fermi 
surface. 



diately see that in the superconducting state, /mnq(f2) 
remains linear in frequency only for 3> A (where 
/mnq(ri) = 7^^). At smaller frequencies, Imllq{^}) van- 
ishes below a threshold at Uth = jA+j + |A_|, where, 
we remind, A± = A (km ± q/2), and km + q/2 and 
km ~ q/2 are discrete pairs of momenta (specified by 
to), which are simultaneously on the Fermi surface. By 
Kramers- Kronig relation, the vanishing of Imllc^{^l) be- 
low the threshold generates a nonzero Rellq{n), which 
comes from frequencies of order A and is therefore part 
of the low-energy theory. This Rellq{il.) dominates the 
spin dynamics below the threshold. 

The structure of nq(ri) in a d— wave superconductor 
has been previously analyzed for q = qo = (tt, tt)—. For 
q = qo, different regions specified by to are all equiva- 
lent, A(kF + qo) = -A(kF) = A, i.e., = 2A. At 
the threshold frequency, ImTlq^^fl) is discontinuous and 
jumps from zero to 7rA7. By the KK relation, Rellqg{^l) 
then diverges logarithmically at fltf^. This divergence 
ensures that for arbitrary couping, Xs(q, ft) — (xj^(ci) — 
nq(ri))~^ has an excitonic pole at some ujres < 2A, where 



nqo (t^r 



Renvoi 



At weak coupling. 



oJres is exponentially close to 2A, and the resonance is 
easily washed out by e.g., disorder. At strong coupling, 
however, the pole is located at small frequencies and is 
weakly affected by disorder. Furthermore, expanding Q 
in powers of fJ, one can easily find that at the lowest fre- 
quencies, nqo(ri) C!C fi^/A, i.e., at strong coupling, when 
i^res *C 2 A, the low- frequency spin susceptibility has a 
magnon-hke form Xs(qo, ^^) « (w^es ~ il^)^^. 

The resonance behavior of nqo(il) sets the crossover 
in the dispersion of an antinodal fermion, for which the 
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scattering by qo is a low-energy process. For a nodal 
fermion, however, the scattering by qo is ineffective, and 
one should analyze other q~^. For a general q 7^ qo con- 
necting two Fermi surface points, we find from Eq. Q 
that the magnitude of the jump in Imllq{^) at the 
threshold frequency Qth — |A(kF)| + |A(kF + q)| is 



S[lmllq{^th 
(1 - sgn{A{kF)A{kF 



^V|A(k^)A(k^ 



-q))) 



(11) 



It then follows that for scattering from a nodal Fermi sur- 
face point ki?^„ to some other point k^^ — kp.n + q along 
the Fermi surface, the jump in ImIlq{Qth) disappears 
because A(kF,n) = 0, even though A(kF,n -I- q) ^ 0. In 
the absence of a jump in Imllq{fl) at flth, Rellq{^l) does 
not diverge when approaching flth from below. Indeed, 
by the KK relation 



dximllq (x) 



Near i7 = flth, this reduces to 

T^^th Jo y + [^th ~ ") 



(12) 



(13) 



When Imllq{il.th + 0+) has a non-zero value (which is 
the case when Imllq{fl) jumps at the threshold), it can 
be pulled out from the integral over y, and Rellq{ilth) 
diverges logarithmically. Without the jump, Imllq{y + 
V,th) vanishes at y = 0+, and the integral over y in l(T^ 
does not diverge. We find from @ that for scattering 
that involves a nodal fermion, ImJlq{y + ilth) oc y^^^, 
hence for < we have from (|13|l . 



Renq{Q) = -7 



(14) 



The minus sign in front of the square-root implies that 
ReHqifl) is negative, i.e., x~^{<lM) = XjHq) ^ ^^i^) 
does not change sign below flth , and the resonance mode 
does not emerge. We recall that a constant, A indepen- 
dent term, has been already pulled out from ReUq{Q), 
hence the negative value is with respect to the normal 
state (as xr^(q) > in a paramagnet, x7^(q) ^ nq(ri) 
is positive for all fl < ilth)- 

Eq. (|14ll can be easily extended to the full complex 
nq(ri) which takes the form 



^2 



(15) 



Note that the square-root functional form of nq(il) sur- 
vives even when the scattering is between a nodal fermion 
and a point not exactly on the Fermi surface. Indeed, 
as long as q is not directed along the zone diagonal, 
one can set, without loss of generality, the velocity of 
a nodal fermion to be along the y axis, i.e., ek = y, 
and the velocity of a fermion at k -I- q to be along x: 
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FIG. 3: a) Different pairs of Fermi surface points separated 
by a momentum qi equal to the node-antinode distance. 
Solid points mark the node (N) and antinode (AN). In ad- 
dition, the other node-antinode wavevector (q2) is shown, 
b) the imaginary and real parts of the particle-hole bub- 
ble xo{<li,^)- The spin polarization operator n(q, f2) = 
Xo^^(q, 0) — Xo^ i'i^^)- Note near discontinuities in Imxo 



1.07A and 1.6A. Here, A 



=40 meV, (5=0.5 



meV, and T=0.5 meV. c) Jmxs (qi, ^) for two values of J(q) 
(x7^(q) = Xo^(<liO) ~ "^(l))- Here, J(qo) was chosen so as 
to yield a resonance at 40 meV for q = qo, and J(qi) is 
- J(qo)(cos(gi2,a) -I- cos(gij;a))/2. 



fk+q = Co + 2: where eo = CkF+q- Substituting this ex- 
pansion into nq(il) cx / Gk,w+oGk+q,w we obtain 

dxdy 



nq(ri) (X j diO 



u) + Q. — y - 
eQ + X + u} 



up- — (eo + — A?{x, y) + i5 



(16) 



Elementary analysis shows that the singular contribution 
to nq(f2) comes from uj = — f2, y = 0, and x — — eoj i-e-j 
from the internal momentum range when both fermions 
are back on the Fermi surface. Furthermore, as typi- 
cal y are infinitesimally small, one can neglect the y- 
dependence of the gap, in which case the momentum 
integral is factorized. Integrating over y, then over lj, 
and finally over a;, one immediately recovers ()15(l with 
Q,th — A(— £0,0), that is the gap at the Fermi surface 
point obtained by projecting k -I- q onto the Fermi sur- 
face along the x direction. 

The square-root behavior of nq(r2) is not the full story, 
however, as the same incommensurate q which connects 
a nodal Fermi surface point with some other k^^ may 
also connect other pairs of Fermi surface points for which 
the superconducting gap is non-zero for both points. If 
the signs of the two gaps are opposite, then, according 
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to Hll() . /mnq(ri) still has discontinuities at the corre- 
sponding threshold frequencies fl'^/^, hence Rellq{fl^fj di- 
verges, and xJ^(q) ^ Rellq{fl) crosses zero at some fre- 
quency below Slj^. We checked this possibility for the 
Fermi surface of optimally doped Bi22\'£^. For node- 
antinode scattering, there are two inequivalent q vectors. 
The smaller of the two, which we label q2, has a dy- 
namic response which is small relative to the larger of 
the two, which we label qi. For qi, the node-antinode 
process at Vtth = ^an is the threshold (an = antinode). 
There are, though, two other inequivalent pairs of vectors 
for which k and k + qi are both on the Fermi surface, 
and A(k) and A(k + qi) have opposite signs. For these 
processes, /r7inq(il) is discontinuous at corresponding 
threshold frequencies fi^^ « 1.07A and 1.6A. We illus- 
trate this in Fig. |3| We found, however, that these extra 
processes do not give rise to a resonance in the spin sus- 
ceptibility for two reasons. First, the dynamic response 
at qi is weaker than that at qo, and the divergence in 
-Rellqj^ is further weakened at the lower energy Q!^^ of 
the two since one of the two k vectors is near the node 
in this case. Thus, the inclusion of any damping (due to 
impurities or finite T) is enough to remove the divergence 
altogether. Second, even in an idealized situation with 
zero fermionic damping, the threshold frequencies fij^ 
for both of these extra scattering processes exceed Vlth, 
hence near £7^^, /mllq^ (f2) is already non-zero. As a con- 
sequence, we find no resonance for Imxsi^x^ ^) (though 
there can be a peak associated with the higher energy 

Summarizing, we argue that for qi which connects 
nodal and antinodal Fermi surface points, there is no ac- 
tual resonance in the spin susceptibility. The imaginary 
part of Xs{(\,^) emerges at Vtth = Aa„ as (Tt - Aq„)^/2 
and has extra bumps at higher energies near threshold 
frequencies for additional scattering processes with the 
same q. Alternatively speaking, the excitonic resonance 
in the dynamical spin susceptibility exists for q near qo 
but gradually vanishes as q approaches qi. The bound- 
ary between these two regions is roughly set by the di- 
agonal node-node scattering vector, which has a length 
intermediate between qo and qi. 



2. fermionic self-energy 

We next compute the fermionic self-energy, Eq. ||2J). 
For an antinodal fermion, both k and k -|- qo are near 
the Fermi surface, and the resonance mode has a strong 
effect on the self-energy. Assuming that Xs {^.^ ^) has a 
magnon-hke form Xs(q, Vt) = xo/i^reAl)-^'^)^ adding a 
small damping term iS to for continuity, and neglecting 
the momentum dependence of uJres both for simplicity 
and because experimentally, the resonance dispersion is 
rather flatS, we obtain from (j^J in the superconducting 



state 



-'SC,an 



(a;) 



3*5sXo f dxdn{uj + n + x) 
Att^vf J + - iuj + - iS 
1 



(17) 



The subscript for the self-energy implies SC = supercon- 
ducting state, an = antinode. We also defined x = VFq± 
and used the superconducting Green's function for free 
fermions 



Lj2 + - A2 - el 



(18) 



The integration over x is straightforward, and performing 
it we obtain 



^SC,an 



(C.) 



3.9^X0 f dn 
u + Vt 



(19) 



This integral is singular near lo = —ujo = — (A + uJres)- 
To see this, introduce uj = —ojq + e, and f2 = ^res + y- 
Substituting these expansions into H19|l and restricting 
to only linear terms in y and e, we obtain after simple 
algebra 



2^SC,Qn(e) 



35fxoVA 



dy 



1 



^'K\/2uJresVF J U + i6 y/-e - y + i5 

(20) 

Splitting the integration over y into integrals over posi- 
tive and negative y and evaluating them separately, we 
obtain 



dy 



1 



2Tr 



y + iS ^-e - y + iS - iS 

Substituting this into (|20|l we obtain, 

Sg^^XoVA V2 



-'SC.an 



SUresVF y/e — iS 



(21) 



(22) 



Separating real and imaginary parts of l/^/T^^JS and 
replacing e back hy lu + luq we finally obtain 



(CJ + CJQ + y/{uJi)+L0Y + (52)1/2 



where 



A 



3g?xo A 



8wfA2 uj^ 



(23) 



(24) 



Note that Ag/ / is dimensionless (with the above definition 
of Xo)- 

We see from (|23ll that the real part of the self-energy 
has a near one-sided singularity. It almost diverges as 
l/^/wo + ^ as LO approaches — from above {\uj\ < 
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Wo), and then rapidly drops beyond — wo, reducing to 
0(6"^ /{to + Wo)) when > ujq. The relation between 
this self-energy and the dispersion is somewhat compli- 
cated in a superconductor, as the MDC lineshape does 
not have a simple lorentzian form because the Green's 
function has ek = VF{kj_ — kp) both in the numerator 
and the denominator: 



G(k,w) = 



ui 



S(w) + Ek 



(W -SM)^- $2(^11 



(25) 



Here $(fc||,a;) is the pairing vertex. It is related to the 
pairing gap A(fc||,Ci;) by $(fcj|,Ci;) = A(fc|| , w)Z(a;), where 
Z{lj) = 1 — i;(w)/w ^S*i2iS4. For simphcity, we neglect 
the frequency dependence of A(fc||,w), i.e., approximate 
A(fc||,a;) by a frequency independent gap A(fc||). Near 
the antinodal points, the gap is near its maximum, i.e., 
is rather flat as a function of , and can be approximated 
by a constant A. Still, the presence of ek in the numer- 
ator of l|25|) implies that the maximum of ImG(\i, uj) is 
shifted somewhat in uj from where the real part of the the 
denominator in H25|l vanishes (in the BCS limit, this effect 
can be attributed to the /c— dependence of the coherence 
factors). This complication, however, affects the form 
of the dispersion mainly for |w| < A, and is less relevant 
near co = —uo where the self-energy is nearly singular. To 
avoid this complication, we neglect the fc— dependence of 
the numerator of the Green's function, and extract the 
dispersion from 



Re 



(26) 



Substituting $(w) in terms of A, and neglecting ImZ{uj) 
(which vanishes for < loq anyway), we obtain from 



(w - j:'scnni^)) Re 



;2- A2 



vpik^-kp) (27) 



Substituting from , we find that the disper- 

sion develops an S-shape for |a;| < wq, precisely as seen 
in the experiments. We illustrate this in Fig. ^ We 
also recall that a near divergence of ReT,an{'^o) imphes, 
by KK transform, a near discontinuity in /mEon(wo), 
both of which give rise to the experimentally observed 
peak/dip/hump behavior of the ARPES intensity— 

For a nodal fermion, the situation is different. A shift 
by qo moves a nodal Fermi surface point to a point where 
the energy is very large, ~ 0.7eV for optimally doped 
_Bz2212. The (tt, tt) scattering is then ineffective. How- 
ever, for q which we analyzed in the previous subsection, 
a nodal fermion can still scatter along the Fermi surface, 
which gives rise to a much larger self-energy. 

The computation of the self-energy J^sc,ni^) requires 
extra care, as we will have to average over all q which 
connect a nodal point with other points on the Fermi 
surface. Besides, even for the integration near a particu- 
lar q, the dispersion of the d— wave gap is essential as it 
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FIG. 4: The self-energy and the dispersion near an antinodal 
point. The coupling is A(qi) = 2. We also set A = cures, i.e., 
iJo ~ 2A, and use a broadening S — 0.3A. Note the 5*— shape 
dispersion near ljq- 



affects the functional form of ImYjsc,n {^) via the soft- 
ening of the singularity in the fermionic density of states 
at fermionic frequencies near the gap at kj? -1- q. As our 
goal is to demonstrate that the self-energy at the nodal 
point does not have the sharp features of the antinodal 
self-energy, we assume for simplicity that (i) the domi- 
nant contribution to Ssc,n(w) comes from node-antinode 
scattering, because of the presence of the density of states 
singularity associated with the antinode, and (ii) that the 
superconducting gap has a flat dispersion in the antinode 
region. 

Substituting the spin susceptibility with H given by 
(|15|l into Eq. ^ and neglecting momentarily the disper- 
sion in Xs (q) around qi , we obtain after integrating over 
momentum near the Fermi surface {VLth — A) 



-'SC. 



n(w) 



A(qi) 



duj'- 



A2 



1 - z^(w + w')2- A2/w,/(qi) 



(28) 



where a;s/(qi) = (7qiXs(qi))^"^- At the lowest frequen- 
cies, expanding to linear order in oj and collecting the 
prefactor, we obtain 



where Asc(qi) = A(qi) 7/'(A/ws/(qi)), and 

dz 1 



(z2 + 1)3/2 i^^^/^r^ 



(29) 



(30) 



such that il^{x) < ip{0) — 1. This implies that the cou- 
pling constant in the superconducting state is somewhat 
smaller than in the normal state. This is in agreement 
with earlier work^2^ At larger frequencies, ReJ^sc.ni^) 
is continuous and reduces to its nearly flat normal state 
form, Eq. H10|l at \uj\ >> A,ujq. The limiting behavior 
resembles that in the normal state, however the crossover 
in H28|) is not analytic, and the self-energy develops a kink 
at = — 2A. This can be most easily seen by evaluating 
the derivative of the self-energy. Indeed, differentiating 
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with respect to u in (|28|l . we obtain 



duj 



duj' 



y(^7)2~A2y(Z7TT7)2~A2 



(1 - iV(^ + c^')'-A2M/(qi))^ 



Near a; = ±2A, the dominant contribution to the integral 
comes from lj' « —l^sgn[ijj), when the two branch-cut 
singularities merge. In this region, uj'{uj + uj') ~ — A^, 
{u + bj'Y K, A^, and 



SC,n 



(CO) 



duj 

-*A(qi)A2 



duj' 



A- 



:(32) 



Evaluating the integral, we find that at, e.g. uj = — 2A, 
dI'mTi{uj) / duj undergoes a finite jump 



dw 



[/mEsc,«(-2A-e) 



dx 



V^Ve-^ 



= 2 



- /mSsa„(-2A + e)] cx 



(33) 



By the KK relation, dRe'S{uj) / dui diverges logarithmi- 
cally at w = — 2A. This behavior is analogous to the 
behavior of the polarization operator with q = (tt, tt) 
near the threshold frequency ilth — 2 A. Evaluating 
dYi'gQ ^{u)) / du) explicitly, introducing lj ~ — 2A -|- e and 
= A -|- x, and expanding to first order in x and in 
e in the two terms in the denominator, we obtain near 
uj = -2A 



duj 



4ws/(qi) 



Im 



dx 



A y/x + i5\/- 



X + iS 
(34) 

Here A ~ A is the upper cutoff for the linear expansion. 
Splitting the integral into the integrals over positive x 
and over negative x, and evaluating them separately, we 
obtain, to logarithmical accuracy 



duj 



-K log- 



2AI 



(35) 



where K = A(qi)A/(2(jJs/(qi)). Integrating this for- 
mula, we obtain 

I]'sc.,„(c.) = Y.'sc,n{^2A) - K (a, + 2 A) log - 



2A| 

(36) 

Note that, contrary to H23|l . the singularity in Eg(^„(w) 
is two-sided, i.e., is symmetric with respect to w -I- 2 A. 
Substituting this self-energy into the dispersion of a nodal 
fermion oj — i}^) = vpik^ — kp), we find that the kink 
in E' gives rise to a kink in the dispersion at = — 2A, 
but the S'-shape form does not emerge. We plot i?eS(a;) 
and the dispersion in Figs.lSJ^ andjSJj. 



ReE(a))/cogf 



(31) 
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FIG. 5: -ReE and the dispersion along the nodal direction. 
a)-f) are obtained with liqiuj) given by dJ. a)-b) are for 
a flat static spin susceptibility Xs(ci) near qi, c)-h) for the 
Ornstein-Zernike form. c)-d) are obtained within Eliashberg 
theory, e)-f) assuming a vanishing Fermi velocity at qi (anti- 
Eliashberg approximation), g) and h) are obtained using the 
spin polarization operator, Eq. I)47|l . with two extra thresholds 
at 0,[f^. In all cases, the dispersion shows a kink around 2 A, 
and does not have the S'-shape form typical for the antinodal 
dispersion (see Fig. 2J. Dashed lines are extrapolations from 
high frequencies; the extrapolated dispersion at oj = has a 
negative Afc, i.e., k inside the Fermi surface. The coupling 
A(qi) = 2 and A/tj^/(qi) = 2. 



The inclusion of the momentum dispersion of Xs(q) 
around qi affects the prefactor, but the logarithmical 
singularity near — 2A survives. Indeed, this singularity is 
just the consequence of the square-root non-analyticity 
in the spin polarization operator, which is confined to 
bosonic frequencies near A. At these frequencies, Ilq^ [uj) 
is small, and the dynamical spin susceptibility can be 
expanded in powers of nq^(w). The momentum de- 
pendence of Xs(q) then just affects the prefactor of the 
(ti;-t-2A) log(a>-f 2A) term in (|35|) . To verify this, we com- 
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puted the self-energy and the dispersion for the Ornstein- 
Zernike form of Xs (q) • In this situation, 



1 



(37) 



(1 - zy(Z7wF^A^M/(qi))i/2 
Again, at low frequencies 

ReY^scA^) = -Asc(qi)w, (38) 
where now A5c(qi) = A(qi) ■(/j(A/cj5/(qi)), and 
dz 1 



(z2 + 1)3/2 (1+2^7^2^)1/2 



(39) 



Near lo = — 2A, we obtain from ()37|l the same functional 
form as in H3t)|) 



E^c. J^) = S^a„(-2A) -k{u, + 2A) log - 



A 



\u + 2A| 

(40) 

where now K — A(qi)A/(4a;s/(qi)). We plot the self- 
energy and the dispersion for the Ornstein-Zernike static 
susceptibility in Figs. |3; and (Sji- We see that the dis- 
persion does not change much from that for a flat Xs {^) , 
namely there is a kink near 2A, but there is no 5— shape 
dispersion as occurs for an antinodal fermion. 

Finally, to verify that the nodal dispersion is not an 
artifact of our computational procedure, we also com- 
puted Tigfj ^{uj) with the Ornstein-Zernike Xs(q) in a 
different (anti-Eliashberg) computational scheme: we as- 
sumed that the Fermi velocity of an antinodal fermion 
nearly vanishes, and integrated over momenta normal to 
the Fermi surface in the bosonic rather than the fermionic 
propagator. In this situation, 

-^AA(ql) , CO' 

log^:££Mz.H^^^2 (41) 
LOsf{q.i) - - <^'Y - A2 

Again, at low frequencies, the self-energy is linear, 

i?eSsc,«(w) = -Asc(qi)w, (42) 

where now Age = A(qi)'0*(A/ws/(qi)), and ^*{x « 
1) = a;|loga;|, and '^*(x » 1) = 7r/4. Near uj = ±2A, 
the derivative dT,sc,n{'-o)/dLu is again singular, but now 



dT,sc,n{^) 
duj 



-AA(qi) / dco' — 



{{uj')^ - + iS) 
1 



^{uj + u;')^~ A^ + iS w,/(qi) - i^{u; + uj')^-A^ 



Expanding, as before near u = —2A and introducing 
UJ ~ — 2A + £ and w' = A -I- x, we obtain from H43|l 



A2A(qi 



dx 



1 



duj 2\/2Atj5/(qi) J-oo x + iS V-e - x + i6 

(44) 

Evaluating the integral using l|21(l we obtain 



d^SC,ni^) 

du 

i.e, near u = — 2A, 



^A2A(qi) 1 
/2Aw^/(qi) + 2A 



(45) 



l^'sani^) = S^c.„(-2A) - K*Rey/2A{iu + 2A) (46) 

where K* — 7rA(qi) A/(u;s/(qi)). As a result, the kink 
at — 2A survives, and E5C',„(— 2A) still does not diverge, 
i.e., there is no S"— shape, antinodal-type dispersion. At 
the same time, the functional form of the non-analytic 
piece changes - the a;log|x| singularity in the real part 
of the self-energy in and (001 j where x = u;-|-2A, gets 
replaced by a one-sided y/x singularity. We plot the self- 
energy, Eq. (|41|l . and the resulting dispersion in Figs.|3^ 
and[Sf. As the two expressions for the self-energy for the 
Ornstein-Zernike Xs(q), Eqs. H37(l and H41|l . represent two 
extremes of large and small fermionic dispersion com- 
pared to the bosonic dispersion, the actual self-energy 
should be somewhere in between, i.e., it is stronger than 
in for \uj\ < 2A, and weaker for \uj\ > 2A. StiU, we 
emphasize that both computational schemes yield a kink 
in the dispersion, but no divergence of i?eS(— 2A) and 
no 5— shape dispersion. 

The actual behavior of the dispersion is more involved, 
as one has to average over all q for scattering from a node 
to some other Fermi surface point. This obviously weak- 
ens the 2A singularity roughly in the same way as the 
singularity in the density of states at A is weakened by 
the momentum dependence of the gap in a d— wave su- 
perconductor. Furthermore, even if by geometrical rea- 
sons, node-antinode scattering at qi dominates the nodal 
self-energy, the actual form of Tiscni^) is more complex 
than in 128|) . I|37|l or H41|) . As we already discussed in 
Sec. Ill B II the square- root behavior of nq(il) persists 
only in a small region near the threshold, while at larger 
f2, additional features in nqj^(ri), associated with the ex- 
istence of two extra pairs of Fermi surface points sepa- 
rated by qi become relevant (see Fig. . 

To qualitatively estimate the relevance of this effect, 
we model Ilqj^ (fi) in Fig.|31by a combination of a square- 
root behavior above f2 = A and a near discontinuity at 
n = 1.6A: 



n" V^^ - A2 + ^(1 + tanh 



(1.6A)2 



(43) 



2^' ■ b 

(47) 

For consistency with Fig. |21 a ^ 1.2 A, and we set 
6'^0.1A(6 = would correspond to a true discontinu- 
ity). In Figs.|3^ andinji, we plot the self-energy and the 
dispersion for this form of nq(il) and flat Xsil)- We see 
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FIG. 6: a) Re E and b) nodal dispersion generated from the 
model of Ref. [TH Circles are the experimental data of Ref. 0| 
for optimal doped Bi2212 at T=40K. The S-shape is replaced 
by a kink due to (1) dominance of the continuum contribution 
to the self-energy relative to that of the resonance and (2) a 
strong energy broadening of 40 meV which was assumed when 
deriving the self-energy. For these calculations, a'res=40 meV, 
o)=40 meV, and the continuum has a gap of 65 meV and a 
cut-off of 500 meV. For this direction, the bare Fermi velocity 
is 3.37 eVA. 



that the dispersion does not change much compared to 
that with just a square-root form of Ilq^ (il) (see Fig.|3)). 
There is still a cusp near 2A and a non S*— shape disper- 
sion. The only new effect is the extension of the crossover 
region above the kink. This is indeed expected as the new 
term in H47|l affects the polarization operator at higher 
frequencies. 

To summarize this subsection, we see that the result of 
various computational procedures is virtually the same: 
the self-energy of a nodal fermion is much less affected 
by superconductivity than the self-energy of an antin- 
odal fermion. The nodal self-energy roughly displays the 
same crossover as in the normal state, from a linear in 
frequency behavior at small frequencies, to a more flat 
behavior at higher frequencies. Superconductivity only 
sharpens the crossover near 2A, but does not give rise to 
any S*— shape features in the dispersion. 



C. a comparison with Ref. Iisl 

It is instructive to compare our results for nodal 
fermions with the earlier study by Eschrig and Normani^. 
These authors also argued that the dominant scattering 
process for a nodal fermion is node-antinode scattering 
by qi (though they included a second process due to 
node-(7r,0) scattering). They used a phenomenological 
form of the spin susceptibility at qi with the resonance 
piece taken as a Lorenztian of width 2a about qo. The 
resulting value at qi is about 8% of that at qg. Added 
to this is a gapped continuum (with a gap equal to the 
threshold value of Imllqg) modeled as a step function 
with an ultraviolet cut-off. This form of Xs(qi, ^) is not 
exactly the one we used above (in our analysis, the reso- 
nance piece is completely absent at qi), but nevertheless 
is rather similar in the sense that a large part of the 
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FIG. 7: a) ReT, and b) antinodal and c) nodal dispersion 
for the electron-phonon interaction. We used Eq. I|49|l for 
the self-energy with uth ~ 2A and coupling constant — 2. 
Solid and dashed lines are for ^/ujth ~ 0.6 and 2, respectively. 
For moderate damping, the dispersion along both directions 
has an S— shape form. Note that the extrapolated dispersion 
crosses the vertical axis at a positive Ak, i.e., for k outside the 
Fermi surface. For the spin-fermion interaction, this crossing 
always occurs at negative Ak for the nodal direction. 



magnetic excitation spectrum is incoherent. Not surpris- 
ingly, the two forms for Xs(qi,'^) yield similar results 
for the dispersion of a nodal electron. In Fig. |H| we plot 
the self-energy and nodal dispersion obtained with their 
susceptibility. We see that this dispersion has a kink, 
but no 5— shape. This was achieved by putting in a sig- 
nificantly large amount of damping which acts to smear 
out the delta function singularity in energy of the reso- 
nance. An S'-shape would still be present, though, if the 
continuum piece is ignored. This result shows that even 
when the resonance is present, the S'-shape of the dis- 
persion emerges only when the resonance contribution to 
S'(w) overshadows the one from the incoherent piece in 
Xs(q, ri). This implies, in particular, that the 5— shape 
dispersion does not emerge immediately away from the 
nodal direction, i.e., it appears somewhere between the 
nodal and antinodal points. 



III. PHONON SCATTERING 

We now perform the same analysis for phonon scat- 
tering. Consider a system of electrons interacting with 
an optical phonon with a frequency ujph and momentum 
q = ^ph ■ The optical phonon propagator can be reason- 
ably approximated by 

X,.(q,c.)^ J^^] (48) 
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where (7(q) rapidly decays at deviations from ciph, and 7 
is the phonon damping rate. The fermionic self-energy 
due to electron-phonon scattering is given by only 
is replaced by g^^. To simplify the discussion, we 
consider a nodal fermion and assume that the Fermi ve- 
locity at ki? + <lph can be neglected. Substituting 
into the self-energy, we obtain 



Re 



Trf3u! 



th 



(49) 



where (3 = g^^ J giq)(Pq and ujth'^ ^ [uoph + A)^ + 
Solving now for the dispersion, we find that at low en- 
ergies, the dispersion is linear with v*p = + A^^) 
and A^^ — n(3/{ujpiiU!^i^). Near uJth, the real part of the 
self-energy for small 7 nearly diverges as lo approaches 
uJth, giving rise to an S-shape dispersion. We illustrate 
this in Fig. ??. The S'-shape could in principle be elim- 
inated if 7 is very large, i.e, the bosonic mode is almost 
overdamped. However, one can easily make sure that to 
avoid the S'-shape, one requires 



4- > I 

th ^ 



ijJ. 



(50) 



This can be relatively easily achieved at weak coupling, 
but for AgP > 1, which is required to fit the low en- 
ergy renormalizaton, the S'-shape is eliminated only at 
nonphysical 7 >> LOth- Furthermore, for Af^ > 2, the 
S'-shape cannot be eliminated for any 7. If we do not 
neglect the Fermi velocity at ki? -I- <\ph , then at the low- 
est frequencies we obtain qualitatively the same result 
as (|23|) . i.e, a nearly one-sided square-root singularity. 
Again, it is very difhcult to get rid of the S— shape form 
of the dispersion at strong coupling without requiring 
that the bosonic mode is totally overdamped. It is there- 
fore very unlikely that the interaction with an optical 
phonon can simultaneously account for a strong Fermi 
velocity renormalization and give rise to a non S— shape 
form of the dispersion. 



IV. CONCLUSIONS AND EXPERIMENTAL 
COMPARISONS 



We conclude therefore that the spin fluctuation sce- 
nario more likely explains the observed features in the 
electron dispersion than the phonon scenario. Within the 
spin fluctuation scenario: (i) in the normal state, there 
is a crossover from a linear to a more steep dispersion at 
around uj^^ , (ii) below the antinodal dispersion devel- 
ops an S— shape form due to interaction with the spin 
resonance with a characteristic energy of A -I- cj^es , (iu) 
the nodal dispersion below Tc develops a kink at 2A, but 
there is no S— shape dispersion as there is no spin res- 
onance for momenta which connect a nodal point with 
other points on the Fermi surface. 



A. a comparison with experiments 

Qualitatively, our picture of an S— shape dispersion in 
the antinodal region below Tc, and a kink dispersion in 
the nodal region which is similar above and below Tc, 
agrees well with the data. Quantitatively, cij<;y(qo) and 
^s}(s\x) relevant for the antinodal and nodal dispersion, 
respectively, were estimated to be ajs/(qo) ~ 20meV 
and Ws/(qi) ~ 40 — 50meV in Bz22132i. Experimen- 
tally, this scale has been detected for the nodal disper- 
sion and is around bOmeV. The resonance frequency uJres 
and the gap A in optimally doped Bi2212 are both near 
40nieV^^, hence the termination of the S— shape in the 
antinodal dispersion and the kink in the nodal disper- 
sion both occur near SOmeV. Experimentally, this scale 
is SOmeV along the antinodal direction and 50 — 70meV 
along the nodal direction^. Both of these gap-related 
scales are smaller in LSCO, but for that material, there 
are no noticeable differences between the normal and su- 
perconducting state dispersions. This implies that the 
effect of the superconductivity on the dispersion is very 
small (though it should be noted that most of the "nor- 
mal state" data for LSCO were actually taken in the 
pseudogap state). We note in passing that the same 
smallness was cited as a reason for the non-observation 
of the resonance peak in LSC02^. 



B. the doping dependence 

Finally, we discuss the doping dependence of the dis- 
persion. There are few systematic studies of the dop- 
ing dependence along the antinodal direction2ii2i2&. Our 
study shows that the magnitude of the S— shape dis- 
persion should increase with underdoping as Xs(qo) in- 
creases, and the coupling constant gets larger. Along 
the nodal direction, the low energy Fermi velocity Vp ~ 
1.6eT^A is relatively doping independent''. In constrast, 
the slope of the high energy dispersion monotonically in- 
creases with underdopingi. Within our theory, the low 
energy Fermi velocity is given by Vp — vp/{l + A(qi)), 
and the coupling constant A(qi) depends on the spin- 
fermion coupling gs and Xs(qi). The coupling constant 
is weakly doping dependent. The susceptibility does de- 
pend on doping via the magnetic correlation length ^, 
but this dependence is non-singular for the node-antinode 
qi 7^ qo- Once |qi — qo| > , the doping dependence 
disappears, and the nodal coupling A(qi) saturates at 
a fixed value. Previous studies by both us and others 
yielded A(qi) ^ 1. This yields a bare velocity ~ 3eVA 
consistent with band theory. 

We believe that the increase of the high energy slope 
with underdoping is a separate effect associated with the 
fact that at high energies, the system progressively de- 
velops SDW precursors. Indeed, at high energies, the 
diagonal scattering by the resonance mode at qo cannot 
be neglected. The argument is simple ~ at high energies, 
the Green's function of an intermediate fermion can be 
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pulled out from the momentum and frequency integral in 
Eq. and the fermionic self-energy acquires the same 
functional form as in the SDW ordered state: 

J^SDwiKu;)^ ^^^^ (51) 

^ - fk+qo 

where Asdw ^ 9s I d^QdflXsil,^) increases with un- 
derdoping. This form is valid for |cl;| ^ loq. Substituting 
this self-energy into the dispersion relation, we find after 
simple algebra that for the negative energies probed in 
ARPES measurements, the maximum of the MDC dis- 
persion shifts from oj « et to cj « et — '-^'sdw/ (^k+qo ~ 
Ek)- As Ek < and ek+qo > 0, the correction to the 
velocity is positive, i.e., the diagonal scattering enhances 
the value of the velocity compared to vp. As Ag^j^y in- 
creases, this effect gives rise to a progressively steeper 
dispersion, in agreement with the data. 

A simple way to appreciate the doping dependence of 



the high energy dispersion is to realize that the MDC 
dispersion for a gapped state is almost vertical for \uj\ < 
^SDW, with some weak dispersion due to coherence 
factorsii. Thus, as the doping is decreased, and a Mott- 
Hubbard pseudogap begins to develop, the high energy 
dispersion is expected to become increasingly steep. 
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